Noise-induced looping on the Bloch sphere: Oscillatory effects in dephasing of qubits 

subject to broad-spectrum noise 
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For many implementations of quantum computing, 1/f and other types of broad-spectrum noise 
are an important source of decoherence. An important step forward would be the ability to back 
out the characteristics of this noise from qubit measurements and to see if it leads to new physical 
effects. For certain types of qubits, the working point of the qubit can be varied. Using a new 
mathematical method that is suited to treat all working points, we present theoretical results that 
show how this degree of freedom can be used to extract noise parameters and to predict a new 
effect: noise-induced looping on the Bloch sphere. We analyze data on superconducting qubits to 
show that they are very near the parameter regime where this looping should be observed. 



PACS numbers: 85.25.Cp, 03.65.Yz, 75.10.Jm 



Motivated by the prospect of quantum computation 
and communication, coherent quantum operation and 
control of small systems has become a central area of 
physics research. The isolation of these systems from ex- 
ternal noise is a key problem, as noise produces decoher- 
ence. In solid-state systems, some level of 1/f or other 
broad-spectrum noise (BSN) is almost always present, 
and is typically difficult to eliminate [l|. Indeed, in su- 
perconducting qubits, single-electron and other tunnel- 
ing devices, this type of noise is recognized as the factor 
chiefly responsible for dephasing 0, y, 0, @| ■ 

One interesting question is the extent to which the 
characteristics of the BSN can be determined by mea- 
surements on the qubit itself. This has been consid- 
ered by several authors @, 0, S, 0- For the most part, 
these authors considered the case of pure dephasing noise. 
Some theoretical work has been done for "mixed" noise, 
which causes both relaxation and dephasing, but this has 
usually been limited to Gaussian approximations flpj ] . 
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asymptotic analysis, or small numbers of RTNs 

In this Letter, we show how to treat mixed noise ana- 
lytically for all times, fully taking into account the non- 
Gaussian effects. This will enable us to show how to back 
out the characteristics of the noise from qubit measure- 
ments. A new physical effect is predicted: noise-induced 
looping on the Bloch sphere. We shall analyze data on 
superconducting flux qubits to show that they are close 
to the regime in which this effect comes into play. How- 
ever, we stress that this effect can occur in any two-level 
system that is subject to BSN, which can include qubit 
implementation from atomic and molecular physics as 
well as solid-state ones. 

The effective Hamiltonian of a qubit is often written 
as H = — i (ecr' z + A<r' x ) — (t) • <?', where e and A 
are the energy difference and tunneling splitting between 
the two physical states. For example, in a flux qubit e 
is proportional to the applied flux through the supercon- 
ducting loop and A is the Joscphson coupling, h' (t) is 
the noise, a random function. We shall be interested in 



the case where hi (t) comes from K random telegraph 
noise sources (RTNs) with a wide range of switching 
frequencies, giving rise to BSN. The coordinate system 
will be rotated by the angle 9 = tan _1 (A/e) so that 
a' z = cos 9a z — sin^cr^ and the qubit energy eigenstates 
are along the z-axis. 9 is called the working point of the 
qubit and H is then 
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where B = y 'e' 2 + A 2 , is the coupling of the fc-th 
RTN to the qubit. Sk (t) switches randomly between the 
values -1 and 1 with an average switching rate 7^. To- 
gether, these sources generate the random field h (t) = 
12k=i Sk (t) 9k- We assume throughout that B S> h{t). 
gk is at an angle 9k to the z-axis, which can be thought of 
as the angle between the noise axis and the energy axis. 
The pure dephasing case, the basis of much work in this 
field, corresponds to 6k = 0. We deal only with classical 
noise, which is generally thought to be appropriate for 
the systems of interest here. Eq[T] is general enough to 
describe almost any qubit subject to classical noise, since 
an arbitrary power spectrum can be constructed with ap- 
propriate choices of coupling constants gk and rates 7&. 
For simplicity, we drop throughout any azimuthal depen- 
dence of the noise and g~k = <7fc(sin#fc, 0, cos 9k). 

In many cases, the RTN sources have a fixed direction. 
For example, it has been shown that for superconducting 
flux qubits the chief noise source is flux noise [l2| , which 
would put the noise along the z-axis (e direction), and 
gives 9u = 9. We also stress that the working point 9 is 
tunable for many different qubit architectures, either by 
changing e or A 13|, 14, 1]|. 

A convenient solution method for noise from an ensem- 
ble of RTNs with an arbitrary distribution of gk and ~/k 
has recently been given [l6[. This " quasi-Hamiltonian" 
method lends itself to perturbative treatment for both 
fast (weak-coupling) RTNs: gk cos 9k <C jk and slow 
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(strong-coupling) RTNs: gkCos9k 7fe. The occurrence 
of the trigonometric factor in these inequalities is crucial: 
it comes from non-analyticity at the point gk cos Ok = "fk 
which leads to qualitatively different behavior in the two 
regimes. The method can treat experiments that involve 
an arbitrary sequence of control pulses. Because it gives 
analytic results for the non-Gaussian effects of BSN over 
the whole range of 9 , we can identify new systematic 
effects. 

In this work we shall restrict attention to two common 
experimental protocols: Energy Relaxation (ER), which 
measures tier (t) = (cr z (t)), and Spin Echo (SE), which 
measures use (t) = (cr x (t)) with a 7r-pulsc at t/2. Addi- 
tional information about the noise may be available by 
the use of more complicated pulse sequences 

The solutions so obtained arc: 
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where M is the total number of slow RTNs and TV is 
the total number of fast RTNs, labeled by m and n re- 
spectively; M + N = K. Ti and are the energy 
relaxation and dephasing rates caused by fast RTNs, 
T3 = Y] m 7m is the overall dephasing rate from slow 
RTNs. e lm = j m / (g m cos9 m ) and e 2m = g m /B Q are 
the small parameters of the perturbation theory. For 
any value of 9, the portion of the spectrum where the 
perturbation theory breaks down is a small fraction of 
the whole, so we expect the approximation to be well- 
controlled. 

The slow RTNs' contribution to the ER signal is much 
smaller than Ti and thus can be neglected. RTN is non- 
Gaussian noise, which accounts for the relatively com- 
plicated appearance of the formulas. We wish to use 
observations of these signals to back out information on 
the distribution of the g's and 7's. 

The fast RTNs' contribution to qubit decoherence is 
fully captured by two decay constants Ti and IV For 
them, our results coincide with Redficld theory [17| and 
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These fast RTNs give rise to purely exponential decay, as 
is well known jl8j |. 

In contrast, the slow RTNs give rise to the r 3 term 
in the exponent of use, which gives exponential decay, 
but also to the more complex and # m -dependent fac- 
tors. These factors give rise to decays which are quasi- 



Gaussian over short range of time and possibly oscilla- 
tions at longer times. 

For purposes of fitting, it is necessary for us to specify a 
not completely general but yet still flexible model for the 
noise. Let n (7) = 2fc=i $ (l ~ 7fc) De the distribution of 
rates and take contant values gk — .<?, 9k = 9. If there is 
a range of couplings then g in the following formulas can 
be regarded as the root-mean-square coupling. We will 
assume a broad noise spectrum by taking n (7) = 017 s-1 
for 7 m j n < 7 < 7 max and otherwise. Note s = gives 
1/f noise. This power-law assumption is often useful for 
analyzing experimental data while the method itself is 
capable of treating arbitrary distributions. The task of 
data analysis is then to determine g,a,s, 7 m j n , and 7 max 
from observations of use (9,t). The data are relatively 
insensitive to 7 ma xj making it difficult to determine. We 
will comment on other limitations of the model below. 
The results of the fitting can be used to make predictions 
for future experiments. 

With these assumptions, the previous results can be 
simplified into 
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where j c = g cos 9 is the critical coupling strength that 
separates fast and slow RTNs and T 3 = J^ c ^,(7)7^7. 
The key point is that by tuning the working point 9, one 
effectively changes the relative number of fast and slow 
RTNs in the environment. Since the two kinds of RTNs 
have different qualitative effects on the qubit, we can use 
this fact to get information about the noise. 
Given these results, it is convenient to define 
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where $5^ (i, 9 = 0) correspond to the "phase-memory 
functional" defined by other authors [Tl|. Note ^2(9) ~ 
Ti(0 = tt/2)/2, thus it drops out in EqO 

We shall now use these formulas to determine some 
parameters of the noise experienced by the flux qubit 
(sample A) in Rcf.[14|. For this system, A/h = 5.445 
GHz and e/h varies from to 1 GHz, meaning that cos# 
varies from to about 0.2. FigQ] shows the fit of Eq[8] 
to the experimental results at cos 9 = 0.18. This gives 
T 3 = 0.99 MHz and y c = 2.1 MHz. The fit is certainly 
excellent, but we note that a Gaussian, as used by the 
authors of Ref.[l4[ fits equally well over this time range, 
with only a hint of deviation at the longest times. The 
working point can be calculated from cos 6* = e/Bq, thus 
fitting 7c(#) determines g/h = 9.6 MHz, as shown in 
FigEJa). Fitting the data of T 3 (9) in FigH;b) gives s = 
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FIG. 1: (Color Online) Echo phase memory functional $se 
data in Fig.4A of Ref . jl3] . We fit the 36 data points to Eq|8] 
(green solid line) and Gaussian model $f £ = e~ Fat (blue 
dashed line) respectively. The center of the first plateau is 
thus predicted to be at rp = 5tt/2"/ c ~ 3.7/zs. 

(true 1/f noise), a = 0.77 and 7 min = 0.048 MHz. Note 
this dependence is very sensitive to the value of s thus 
can be used to identify the noise distrubition n("f). We 
have carried out the analysis for the results of a similar 
experiment of Ref. [l5| . All results are given in Table [U 
An interesting consistency check is to compute Ti from 
the model with the parameters determined since Ti can 
be fitted from nERif) independently. For the data of 
Ref. [3 (column A of Tabled, we find T ( * h) = 0.02 MHz, 

while the measured value is Tj 6 = 0.65 MHz. This 
strongly suggests an additional source of high-frequency 
noise not captured by our noise model. For the data 
of Ref.0 (column B), we find ri" l) = 5.4 MHz and 
T^ x ' = 5.7 MHz; this shows that the noise model is 
reasonably complete for this experiment. 

The experiments are typically fit by the Gaussian the- 
ory for dephasing noise 0. Galperin et al. 
have pointed out the importance of non-Gaussian effects 
and have shown numerically that these effects appear in 
&SE (t) at strong coupling and longer times. This takes 
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A/h (GHz) 


5.445 


3.9 


e/h (GHz) 


- 1 


- 0.925 


T[ ex) (MHz) 


0.65 


5.7 


g/h (MHz) 


9.6 


135 


j min (MHz) 


0.048 


0.098 




0.77/7 


0.75/7 


Yf l) (MHz) 


0.02 


5.4 



TABLE I: Noise characteristics extracted from Ref. [3 (col- 
umn A) and Ref.[lB] (column B). 




FIG. 2: Fitting of sample A data from Ref. [3. Both 7 C and 
T3 are fitted from &se at various working point, (a) Critical 
rate y c versus the working point cos 6. (b) Linear regression 

to r 3 = Q (7 C — Jmin)- 

the form of "plateaus" . Eqs [7] and [8] allow us to quantify 
and systematically analyze these deviations from Gaus- 
sian behavior at arbitrary 9. 

Note that for BSN, a is the parameter that controls 
the appearance of the plateaus for 1/f noise. Physi- 
cally, it describes the overall scale of the noise since 
K = J lmaa: n(j)d r y. Fig[3] shows the behavior of $se (t) 
for various values of a. In these plots we have used the 
parameters of Ref. [3 and varied only a. The experi- 
mental value a = 0.77 is closest to FigfSJb). Fig[3{c) 
shows that the experiments are not far from seeing some 
non-Gaussian behavior. 

FigfS] shows that as a is reduced a series of plateaus 
becomes visible and the phenomenon is better described 
as oscillations, as would be expected from EqJS] This 
could be done by reducing the number of RTNs, which 
is presently a very active area of research [2(| • FigfSfd) 
shows these plateaus quite clearly, and confirms that the 
period of these oscillations is given by 2ir/(g cos 6). The 
observation of a period that is adjustable by changing 
the working point would be an unambiguous signature of 
the effect. 

The physical interpretation of the oscillations is sim- 
ple. $se if) is a measurement of the probability that 
the Bloch vector has returned to its starting point on 
the equator in the rotating frame. The Gaussian theory 
treats the sphere as a plane and the decoherence as a 
random walk on that plane. In the more rigorous theory, 
there is an additional process - a return to the starting 
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FIG. 3: (Color Online) <3?se with different noise intensity a. 
In all subfigures, solid black line corresponds to cos# = 0.2 
and solid red line corresponds to cos (9 = 0.5. The correspond- 
ing calculated from Gaussian approximation are plotted 
as dashed green lines. Note Fq oc ^/acos8 and Gaussian 
approximation works well in the short time region until the 
emergence of the first plateau, if there is one. All subfigures 
share the same axes labels as in (c). 

point after looping the sphere. This process is only pos- 
sible for the slow RTNs whose fields can persist over a 
time long enough for the loop to occur. It is due to mo- 
tion along lines of latitude on the sphere and only the 
z-componcnt of these fields produces this. This explains 
why the characteristic looping time is 27r/(c/ cos (9). 

This looping depends on having not too much spread 
in the distribution of the g*k- Clearly a very wide spread 
would wash out the oscillations, as seen in Eqj3] On the 
other hand, a wide spread that weights the frequency 
components of the power spectrum unequally is not con- 
sistent with the usual picture of the origin of 1/f noise as 
being due essentially to a spread in 7^. 

In conclusion, we have exploited a new mathematical 
method to obtain the qubit decoherence behavior as a 
function of working point 9 for qubits that are subject to 
BSN coming from classical sources. Varying 9 changes 
the ratio of strongly and weakly-coupled RTNs. It there- 
fore furnishes a convenient method to back out noise 
parameters from observations of the qubits themselves. 
This allows us to refine the picture of non-Gaussian os- 
cillations in qubit decoherence. In particular we have 
supplied a method to observe the oscillations, and have 
identified their physical origin as noise-induced looping 
on the Bloch sphere. This looping should be observable 
in many different qubit implementations. 
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